The Reactor Antineutrino Anomaly 
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Recently, new reactor antineutrino spectra have been provided for 235 U, 239 Pu, 241 Pu, and 238 U, 
increasing the mean flux by about 3 percent. To a good approximation, this reevaluation applies 
to all reactor neutrino experiments. The synthesis of published experiments at reactor-detector 
distances < 100 m leads to a ratio of observed event rate to predicted rate of 0.976 ±0.024. With our 
new flux evaluation, this ratio shifts to 0.943±0.023, leading to a deviation from unity at 98.6% C.L. 
which we call the reactor antineutrino anomaly. The compatibility of our results with the existence of 
a fourth non-standard neutrino state driving neutrino oscillations at short distances is discussed. The 
combined analysis of reactor data, gallium solar neutrino calibration experiments, and MiniBooNE- 
v data disfavors the no-oscillation hypothesis at 99.8% C.L. The oscillation parameters are such that 
|Am 2 ew | > 1.5 eV 2 (95%) and sin 2 (26> ncw ) = 0.14 ± 0.08 (95%). Constraints on the 0i 3 neutrino 
mixing angle are revised. 
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I. INTRODUCTION 

Neutrino oscillation experiments over the last twenty 
years have established a picture of neutrino mixing and 
masses that explains the results of solar, atmospheric and 
reactor neutrino experiments [ij. These experiments are 



(90% C.L.) for IAto^I ~ 2.4 1CT 3 eV 2 . Forthcoming 
kilometer baseline experiments with multiple detectors 
are underway to provide a clean measurement of the last 
undetermined neutrino mixing angle 0±3 

EMI. 

The KamLAND experiment detected electron antineu- 
trinos from surrounding reactors located at an average 



consistent with the mixing of v e , Vu and v T with three distance of 180 km, thus probing the Am^ - 10 4 eV 



mass eigenstates, v\, V2 and 1/3. In particular, the mass 
differences are required to be lAm^l ~ 2.4 10 -3 eV 2 and 
Amli/IAmfuJ ~ 0.032. 

Reactor experiments have played an important role in 
the establishment of this pattern. In experiments at dis- 
tances < 100 m from the reactor core, at ILL-Grenoblc, 
Goesgen, Rovno, Krasnoyarsk, Savannah River and 
Bugey [ll| , the measured rate of v e was found to be 
in reasonable agreement with that predicted from the re- 
actor antineutrino spectra, though slightly lower than ex- 
pected, with the measured/expected ratio at 0.976±0.024 
(including recent revisions of the neutron mean lifetime). 
This opened the way to longer baseline experiments mea- 
suring the v e survival probability P ee at various dis- 
tances. 

Including the three known active neutrinos, P ee can be 
written in first approximation as 



P£{L) = 1 - sm z {20 i )s\TL Z 1.27 



.Am 2 [eV 2 ; 



E Pe [MeV] 



(1) 



where Ar 



Ami and 0i ~ 9\2 if the baseline exceeds 



13 



a few tens of kilometers, and Am- = Am 2 ^ and 0i = 
if it does not exceed a few kilometers [l2| . 

The CHOOZ experiment was located 1 km away from 
two 4.3 GW reactors 13] , and did not report any neutrino 
oscillation in the parameter region Am| x > 10~ 3 eV 2 . 
In addition CHOOZ data provides the world's best con- 
straint on the 6*i3 mixing angle value, at sin 2 (2&13) < 0.14 
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region. KamLAND successfully reported a strong oscilla- 
tion signal in agreement with solar neutrino data [l7j . A 
three-flavor oscillation treatment shows that KamLAND 
also has a weak sensitivity to the #13 mixing angle, as 
P ee ~ cos 4 0i 3 P^, leading to the anti-correlation of 
sin 2 613 and sin 2 6*12- Recent global fits including solar 
neutrino and reactor oscillation data indicated a prefer- 
ence for non-zero #13 at about 1.5 er [l8l [l9j. 

In preparation for the Double Chooz reactor experi- 
ment, we have re-evaluated the specific reactor antineu- 
trino flux (^/fission), improving the electron to antineu- 
trino data conversion [2(j] ■ The method relies on detailed 
knowledge of the decays of thousands of fission prod- 
ucts, while the previous conversion procedure used a phe- 
nomenological model based on 30 effective beta branches. 
Both methods are constrained by the well-measured ILL 
spectrum of fission induced electrons that accompanies 
the antineutrinos [2"lrl24j ]. 

The new calculation results in a small increase in the 
flux of 3.5%. Although this increase has no signifi- 
cant effect on KamLAND's solar parameter results, when 
combined with the previously reported small deficits at 
nearer distances, it results in a larger average deficit 
of 5.7%, at 0.943 ± 0.023. We call it the reactor an- 
tineutrino anomaly, significant at the level of 98.6% C.L. 
This deficit is consistent with being independent from the 
distance to the reactor core at distances > 15 meters, the 
distance of the Bugey-3 experiment [3]. 

If this deficit is due to neutrino mixing, it could be ex- 
plained by an energy-independent suppression of the v e 
flux at distances > 15 meters. This requires a neutrino 
with a |Am 2 cw | > 1 eV 2 . The mixing amplitude with the 
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v e must be sin 2 (26* ncw ) ~ 0.115. The required |Am 2 ew | 
is significantly larger than those required by solar and 
atmospheric experiments. This suggests, if the neutrino 
mixing hypothesis is the correct explanation, the exis- 
tence of a fourth neutrino, beyond the standard model. 
We would like to stress here that other explanations are 
also possible, such as a correlated artifact in the experi- 
ments, or an erroneous prediction of the antineutrino flux 
from nuclear reactor cores. 

We begin by discussing the use of new antineutrino 
spectra and provide a revised value of the predicted cross 
section per fission, in Section [TTJ We then revisit the 
short baseline oscillation results (< 100 m) and CHOOZ 
in Section IIII| revealing a reactor antineutrino anomaly, 
which is discussed in Section IIVI The compatibility of 
our results with the existence of a fourth non-standard 
neutrino state is discussed in Section [V] The CHOOZ 
and KamLAND sensitivities to the #13 mixing angle are 
revisited in Sections IVI Al and IVI Bl Their combination 
is discussed in Section lVI CI Two experimental programs 
for testing the anomaly are discussed in Section IVlII We 
conclude in Section IVIIII 



II. NEW PREDICTED CROSS SECTION PER 
FISSION 

Fission reactors release about 10 20 v e GW~ s _1 , which 
mainly come from the beta decays of the fission products 
of 235 U, 238 U, 239 Pu, and 241 Pu. The emitted antineu- 
trino spectrum is then given by: 

Stot{E v ) = Y,fkS k {E v ), (2) 

ft 

where fk refers to the contribution of the main fissile nu- 
clei to the total number of fissions of the k th branch, and 
Sk to their corresponding neutrino spectrum per fission. 

For the last 25 years the v e spectra have been esti- 
mated from measurements of the total electron spectra 
associated with the beta decays of all fission products of 
235 U, 239 Pu, and 241 Pu. Thin target foils of these iso- 
topes were irradiated with thermal neutrons at the ILL 
reactor [22T4241 ] . The measured spectra then had to be 
converted from electron to antineutrino spectra invoking 
a set of 30 effective beta-branches, adjusted to reproduce 
the total electron spectrum [39| . 

Recently we revisited the conversion procedure with a 
novel mixed-approach combining the accurate reference 
of the ILL electron spectra with the physical distribution 
of beta branches of all fission products provided by the 
nuclear databases [2C| . This new approach provided a 
better handle on the systematic errors of the conversion. 
Although it did not reduce the final error budget, it led 
to a systematic shift of about 3% in the normalization 
of 235 U, 239 Pu, and 241 Pu antineutrino fluxes, respec- 
tively. This normalization shift has been attributed to 
two main systematic effects in the original conversion of 
the ILL electron data. At low energy [E v < 4 MeV) 



the implementation of Coulomb and weak magnetism 
corrections to the Fermi theory in the new approach 
turned out to deviate from the effective linear correction 
(0.65 x [E v — 4 MeV) in %) used in the previous work. 
At high energy [E v > 4 MeV), the converted antineu- 
trino spectra become very sensitive to the knowledge of 
the charge Z of the nuclei contributing to the total spec- 
trum. In the previous approach, only the mean depen- 
dence of Z versus the end-point of the effective beta- 
branches had been used while in the new conversion we 
had access to the complete distribution, nucleus by nu- 
cleus. These two effects could be numerically studied and 
confirmed on various independent sets of beta-branches. 

Because 238 TJ nuclei undergo fission with fast neutrons, 
the associated electron spectrum could not be measured 
in the thermal neutron flux of the ILL reactor. Therefore 
the ab initio summation of the u e from all possible beta 
decays of fission products was performed to predict the 
neutrino spectrum [251 ] - In Ref. [201 ] we provided a new 
prediction with an estimated relative uncertainty of the 
order of 15% in the 2-8 MeV range. This uncertainty of 
ab initio calculations is still too large to be generalized 
to all isotopes but it is sufficiently accurate in the case 
of 238 U, which contributes to less than 10% of the total 
fission rate for all reactors considered in this work. An 
ongoing measurement at the FRM II reactor in Garching 
will soon provide experimental constraints [26j. 

When predicting the antineutrino spectrum of a re- 
actor Stot(Eu), a weighted sum of the four antineutrino 
spectra Su{E v ) is computed according to the considered 
fuel composition, which can be different in each exper- 
iment. The object of this article is to analyze the im- 
pact of the above-mentioned ~ 3 percent shift on past, 
present, and future experiments. 

Generally reactor neutrino oscillation experiments 
search for the reaction: 

v e + p — > e + + n, (3) 

where an electron antineutrino interacts with a free pro- 
ton in a detector, often filled with scintillator. The re- 
action cross section can be precisely computed with the 
V- A theory of weak interaction [2{| : 

fJV-A(-Ee) = KPeE e (l + S rcc + S wm + (5 rad ), (4) 

p e and E e being the momentum and energy of the 
positron, and the 5 being the energy dependent recoil 
(<5 rcc ), weak magnetism (S wm ), and radiative (£ ra d) cor- 
rections. On the one hand, the prefactor k can be written 

K= G|^^ « 2 

where Gf is the Fermi constant, 9c the Cabibbo angle, 
A^ ne == 0.024 the inner radiative corrections taken from 
Ref. [28|, and A = 1.2694 the form factor ratio of the 
axial to vector coupling constant. Parameters are taken 
by default to their latest PDG average values [l[ . Using 
this parameterization Ref. [301 ] obtained a reference value 
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of k — 0.952 10~ 43 cm 2 , using A = 1.2670 close to the 
average value of the year 1999 (31J. This value must now 
be updated to 1.2694 ill], leading to k = 0.955 10~ 43 cm 2 . 

On the other hand, Eq. 2] can be normalized to the 
/3-decay of the free neutrons and k can be written as: 



2^ 2 



5 fR n 



m 5 J 



(0) 



r n being the neutron mean lifetime, and / 
1.71465(15) the phase-space factor for beta-decay of the 
free neutron taken from Ref. [27|, including outer ra- 
diative corrections. Over the last 15 years the neutron 
mean lifetime has evolved, from 926 s (value used in the 
ILL experiment Q) to its current PDG average value of 
885.7 s [1]. Based on this parameterization the current 
prefactor k is 0.956 10~ 43 cm 2 . This is our default value 
in this publication. For the computation of the differen- 
tial cross section we used [2{|- Our results agree within 
0.1% with the results published in Ref. [3(3] that super- 
sede those of Ref. [3!| . 

It is worth mentioning that there is an ongoing con- 
troversy about the world average neutron lifetime, with 
Ref. [Ij] finding r„ = 878 s, and making the world aver- 
age 'suspect'. A lower value has also been found by the 
MAMBO-2 group [3S| . The new world average should 
then evolve and settle to 881.4(1.4) s in 2011 (Ref. [H 
and private communication from K. Schreckenbach) . 

We note here that the average value of A may depend 
on neutron lifetime measurements. However other exper- 
iments studying angular correlations between the neu- 
tron spin and the emitted electron provide independent 
measurements of A. Recently Ref. [32J finds 1.27590 and 
Ref. H3 favors A = 1.2750(9). A value of A = 1.274 
would be consistent with the latest neutron life time mea- 
surements, to be averaged to r„=881.4 s [35}. Thus the 
latest values of A and r n point to a reevaluation of the 
prefactor k to 0.961 10~ 43 cm 2 . With this change, the 
cross section would increase by 0.5% compared to the 
calculations in this work. 

The outgoing positron and incoming antineutrino en- 
ergies are related by 

A EJE e + A) 1 (A 2 - m 2 ) 

where A = M n — M p [29|, [3t| . The prediction of the cross 
section per fission is defined as: 

/>oo 

^r d= J w£>v- A (£,)^=^/ fe ^ fe ed , (8) 

where 5tot is the model dependent reactor neutrino spec- 
trum for a given average fuel composition (/&) and ov-A 
is the theoretical cross section of reaction ((3]). The <J^ d 
are the predicted cross sections for each fissile isotope. 
Experiments at baselines below 100 m reported either 
the ratios of the measured to predicted cross section per 
fission, or the ratios (R) of the observed event rate to the 
predicted rate. 



Accounting for new reactor antineutrino spectra [2G | 
the normalization of predicted antineutrino rates, erj ^ , 
is shifted by +2.5%, +3.1%, +3.7%, +9.8% for k= 23 ^U, 
239 Pu, 241 Pu, and 238 U respectively (see Table Q}. In 
the case of 238 U the completeness of nuclear databases 
over the years largely ex plai ns the +9.8% shift from the 
reference computations [25J. The new predicted cross 
section for any fuel composition can be computed from 
Eq. ©. By default our new computation takes into ac- 
count the so-called off-equilibrium correction of the an- 
tineutrino fluxes (increase in fluxes caused by the decay 
of long-lived fission products). 



III. IMPACT ON PAST EXPERIMENTAL 
RESULTS 

In the eighties and nineties, experiments were per- 
formed at a few tens of meters from nuclear reactor cores 
at ILL, Goesgen, Rovno, Krasnoyarsk, Bugey (so called 3 
and 4) and Savannah River 0-0, E[ ■ Following these pi- 
oneering results middle- and long-baseline experiments 
were performed at CHOOZ [H and KamLAND ^ in 
the late nineties and after. In this section we revised 
each ratio of observed event rate to predicted rate. The 
observed event rates and their associated errors are un- 
changed. The predicted rates are reevaluated separately 
in each experimental case. 



A. The Bugey-4 integral measurement 

The Bugey-4 experiment [3j performed the most 
precise measurement of the inverse beta decay cross 
section, obtaining tr^ ugGy =5. 752+0. 081 in units of 
10~ 43 cm 2 /fission. This experimental result was com- 
pared to the predicted cross section per fission, <jP rcd ' olc \ 

Using Ref. Q inputs, and the former converted 

neutrino spectra from [22|, [4(J El[ (except for the 
238 U neutrino spectrum taken from 20]), we get 
<r? re '° =5.850 10~ 43 cm 2 /fission, in good agreement 
with Bugey-4's estimation (see Table [I]). We used 
the average fuel composition of Ref. Q, 235 U=53.8%, 
239 Pu=32.8%, Pu 241 =5.6% and U 238 =7.8%, in fractions 
of fissions per isotope. We note here that the pub- 
lished Bugey-4 cross section calculation does not account 
for the contribution of long-lived fission products (off- 
equilibrium effects). The reference electron spectra used 
for the determination of S^ot were acquired after an irra- 
diation time in the ILL reactor flux of less than 1.5 day for 
all isotopes. But in the Bugey experiments at a commer- 
cial PWR the irradiation time scale was typically about 
one year. Stot should thus be corrected for the accu- 
mulation of long-lived fission products in the low energy 
part of the spectrum. We had to turn off these effects in 
our computation to recover the (T P rcd ' old value. Including 
these effects in our calculations would lead to a +1.0% 
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r,lrl r^i 


I16W 


prcd 

°7, 235 u 

prcd 
°/,239Pu 
pred 

cr /, 238 u 

prcd 
°"/,241Pu 


fi 39+1 9% 
4.19±2.4% 
9.21±10% 
5.73±2.1% 


fi CI _|_o 1 1 % 

4.34±2.45% 
10.10±8.15% 
5.97±2.15% 


prcd 

f 


5.824±2.7% 


6.102±2.7% 


a^ UECy /a p / cd 


5.752±1.4% [3J 
0.987±1.4%±2.7% 1 0.943±1.4%±2.7% 



TABLE I. Individual cross sections per fission per fissile iso- 
tope, &Yk d - The slight differences in the ratios are caused by 
differences in off-equilibrium effects. Results of the compar- 
ison of a^ ugcy at Bugey-4 in units of 10 -43 cm 2 /fission with 

the former and present prediction, o"5 red , are shown. 



discrepancy. An over or under-estimation of the irradia- 
tion time by 100 days would lead to a systematic error 
on the off-equilibrium correction below 0.1%. 

Computing the ratio of observed versus expected 
events i?Bu ge y=o'^ ugey /cr^ red , Bugey-4 concluded that 
there was good agreement with the V-A theory of weak 
interactions, and that the neutrino flux emitted by PWR 
reactors is sufficiently understood to be computed using 
the parameters of [22|, 0, Hl[ . 

Applying the new reactor antincutrino spec- 
tra we obtain a new value of the cross section 
per fission of cr^ Icd ' now = 6.102 ± 0.163 in units 
of 10~ 43 cm 2 /fission. We thus revised the ra- 
tio i? Bug oy-4=0.943 ± 0.013 (stat + syst) ± 0.025 (S to t), 
which is now 2.0 standard deviations from the expecta- 
tion without neutrino oscillations. This creates a tension 
between the measurement at Bugey-4 and the new 
predicted cross section per fission. In the next sections 
we will see that the other reactor neutrino experimental 
results confirm this anomaly. 

B. The ILL neutrino experiment 

In the early eighties the ILL experiment measured the 
antineutrino induced positron spectrum at a distance of 
8.76 m from the core of the Laue-Langevin fission reac- 
tor. Its fuel is highly enriched uranium (93% 235 U), sim- 
plifying the computation of the predicted antineutrino 
spectrum rate and shape. In Ref. 2] the ILL experiment 
reported a ratio of measured to predicted event rates of 
■Rill,80=0.955 ± 0.035 (stat) ± 0.11 (syst + S tot ). 

In 1995 this result was revised by part of the ILL 
collaboration [36j . The 1980 reactor electron spectrum 
for 235 U of Ref. [2l| was re-measured in 1984 by 
the same authors as [2lJ. The neutron mean lifetime 
was corrected from 926 s to 889 s, increasing the pre- 
dicted cross section by 4%. Moreover in 1990 it was 
announced that the operating power of the ILL reac- 
tor had been incorrectly reported at the time of the 



neutrino experiment, by +9.5%. This reanalysis led to 
-Rill,95 =0.832 ± 0.029 (stat) ± 0.0738 (syst -I- S tot ), ex- 
cluding the no-oscillation case at 2 a. 

According to the new spectra of Ref. 20] the an- 
tineutrino rate is increased by ^3.5% (see Table In)) . 
A slight neutron mean lifetime correction leads to 
an additional +0.37% shift. The new ratio is thus 
fliLL,new=0.802 ± 0.028 (stat) ± 0.071 (syst + S tot ), sig- 
nificantly deviating from its expected value. 



C. Bugey-3, Goesgen, Krasnoyarsk, Rovno, SRP 

We now study the impact of the new antineutrino spec- 
tra on experiments at Bugey (called Bugey-3), Goesgen, 
Rovno, and Krasnoyarsk |^-|9(, which measured the re- 
actor antineutrino rate at short distances, between 15 m 
and 95 m, but less accurately than Bugey-4 Jj|. 

Accounting for new reactor antineutrino spectra (2pj 
the Bugey-3 [J] ratios of observed versus expected events 
between 1 and 6 MeV decrease by 3.7%. However, we 
have to apply a further correction to account for off- 
equilibrium effects. Assuming 300 days of irradiation this 
leads to an additional increase by +1.0%. The Bugey-3 
ratios become: 0.946 ± 0.004 (stat) ± 0.048 (syst + S to t) 
at 15 m, 0.952 ± 0.01 (stat) ± 0.048 (syst + S to t) at 40 m 
and 0.876 ± 0.126 (stat) ± 0.048 (syst + Stot) at 95 m. 
Note that uncertainties on Stot are included in the er- 
rors quoted by the Bugey-3 collaboration. 

A similar analysis is performed with the Goes- 
gen Q results. The new Goesgen ratios shift 
to 0.966 ± 0.017 (stat) ± 0.060 (syst + S to t) at 38 m, 
0.991 ±0.019 (stat) ± 0.062 (syst + S tot ) at 45 m and 
0.924 ± 0.033 (stat) ± 0.062 (syst + S tot ) at 65 m. 

In the eighties a series of experiments were performed 
at the Rovno nuclear power station in the Soviet Union. 

In 1988 the collaboration published measurements 
at 18 m and 25 m from the reactor core Q. Five 
measurements were performed with two different de- 
tectors: an integral detector to measure the absolute 
cross section per fission of inverse /3-decay (labeled II 
and 21), and a scintillator spectrometer to measure 
both the absolute and differential cross section per fis- 
sion (labeled IS, 2S, and 3S). The neutron lifetime 
was taken to be 898.8 s. We accounted for the av- 
erage fuel composition for each run which was pub- 
lished in [6j. As an indication, the average fuel com- 
position over the five measurements is 235 U=59.6%, 
239 Pu=28.3%, 241 Pu=4.6% and U 238 =7.5%.Thc five re- 
sults of the observed over predicted event ratios are re- 
ported in table [Til At 18 m the average ratio is shifted 
from i?Rovno88,i8m,ncw=0.995 ± 0.060 (stat + syst + Stot) 
to i?Rovno88,i8m,ncw=0.944 ± 0.057 (stat + syst + Stot)- 

In 1991 the Rovno integral experiment 0] published 
a cross section per fission of (J J' ovn ° 91 = 5.85 ± 0.17 in 
units of 10 3 cm 2 /fission, 18 m away from a nuclear 
core with an average fuel composition of 235 U=61.4%, 
239 Pu=27.4%, 241 Pu=3.8% and U 238 =7.4%. They 
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predicted the cross section cr^ r ^ 1 J^ 91 = 5.94 ± 0.16 
in units of 10~ 43 cm 2 /fission, and thus obtained 
the ratio i?R O vno9i,oid=0.985 ± 0.037 (stat + syst + Stot)- 
We recomputed the cross section per fission ac- 
cording to the new antincutrino spectra and found 
CT /Rovno9i = 6-223 ± 0.17. The new ratio is thus revised 
to' i? R ovno9i,now=0.940 ± 0.036 (stat + syst + S tot ). We 
note that the correction of the neutron mean lifetime 
contributes 1.8% to the shift of the ratio. 

In 1984 a neutrino experiment operated at the 
Krasnoyarsk reactors 0], which have an almost pure 
235 U fuel composition leading to an antineutrino spec- 
trum within 1% of pure 235 U, and operate over 
50 day cycles. They measured the cross section per 



ferent averaged fuel composition according to: 



fission at two distances, a 



33 m and a 



Krasno,92m 
/ 



Krasno,33m 

/ 



= 6.19 ±0.36 at 



6.30 ±1.28 at 92 m, in units 
of 10 -43 cm 2 /fission. They compared it to the pre- 
dicted cross section of 6.11 ± 0.21 10~ 43 cm 2 /fission, 
based on the Ref. (21] 235 U measurement instead of 
Ref. [221 ]. Correcting the neutron mean lifetime and 
using the new antineutrino spectra we obtain a pre- 
dicted cross section of 6.61 10 -43 cm 2 /fission, assum- 
ing a pure 235 U spectrum. This leads to the ra- 
tios i?Krasno, 33m=0. 936 ± 0.054 (stat + syst ± S'tot) and 
J?Krasno, 92m=0. 953 ± 0.195 (stat ± syst + ^tot), at 33 m 
and 92 m, respectively. In 1994 two other measure- 
ments were performed 57 m from the Krasnoyarsk re- 
actors They measured cr^ rasno ' 57m = 6.26 ± 0.26 
at 57 m, and compared it to their predicted cross 
section of 6.33 ± 0.19 10~ 43 cm 2 /fission, based on 
Ref. [22j and in agreement with our reevaluation 
using previous reference antineutrino spectra. Us- 
ing the new values of Ref. [2(| we revise the ratio 

i?Krasno,57m=0.947 ± 0.047 (stat + Syst + Stot)- 

From the neutrino pioneering experiments led by 
F. Reines and C. Cowan [loj to the nineties, a series of re- 
actor antineutrino measurements were performed at the 
Savannah River Plant (SRP), a U.S production facility 
for tritium and plutonium. For neutrino energies between 
2 and 8 MeV the spectrum difference between SRP and a 
similar core with pure 235 U fuel was estimated to be less 
than 0.5%. We make use of the latest results published 
in Ref. [01 ■ Measurements were reported at two different 
baselines, 18 m and 24 m. The new SRP ratios are reeval- 
uated to 0.953 ±0.006 (stat) ±0.0353 (syst + S tot ) and 
1.019 ± 0.010 (stat) ± 0.0377 (syst + S tot ), respectively. 



D. CHOOZ and Palo Verde 



Based on the good agreement between cr P lcd > old anc [ 



cr^ ugey obtained at Bugey-4 Q, the CHOOZ experi- 
ment [13( decided to use the total cross section per fission 
measured at Bugey-4, correcting a posteriori for the dif- 



Chooz 



Bugey . V - ^ / f Ch 



7 



+ 



k 



f, 



Bugcy \ prcd,old 



f,k 



, (9) 



where /^ hooz and /^ ugoy are the contributions of the k th 
isotope to the total amount of fissions at the CHOOZ 
and Bugey-4 experiments [37|. This explicitly means 
that the expected number of events was computed us- 
ing er^ ugcy rather than cr P rod > old , thus absorbing a -1.3% 
difference on the overall normalization (see Table IJ). 
This also has the effect of reducing the error on the 
neutrino detection rate from 2.7% to 1.6% (including 
the uncertainty on the fission contributions /&). Ac- 
counting for an uncertainty on off-equilibrium effects, 
CHOOZ quoted a final neutrino spectrum error of 1.9%. 
As shown above, the new values from Ref. [2(1 lead to 



r prcd,new 

an increase of o~ f = 



1.048 a 



prcd.old 

f 



With this 



sizeable discrepancy between measured and computed 
cross sections, the CHOOZ experiment cannot rely any- 
more on the effective cross section per fission measured 
at Bugey-4, assuming no-oscillation at baselines of less 
than a few tens of meters. We thus revise the ratio 
to i? C hooz=0.961 ± 0.027 (stat) ± 0.032 (syst + S tot ) (see 
Table Hlf . where the 3.3% systematic error now includes 
the 2.7% uncertainty on cr P rcd ' ncw (See Appendix |A"|). 
A crude analysis of th e impact of Ref. I20I 
Palo 



on raio Verde 
fication of the 
sus expected v, 

^PaloVerde 



data 
average 
by 



44] 



leads to the modi- 
ratio of detected ver- 
roughly —3.5%, leading to 
0.975 ± 0.023 (stat) ± 0.055 (syst + S tot ). 



We also note that Palo Verde's uncertainty would 
increase from 5.3% in Ref. [44[ to 5.6% according to our 
prescription. 



IV. THE REACTOR ANTINEUTRINO 
ANOMALY 

A. Rate information 

Up to now we independently studied the results of the 
main reactor neutrino experiments using a new value 
of the cross section per fission, cr P rcd > ncw _ The ra _ 
tios of observed event rates to predicted event rates, 
R = Aobs/A pro d, are summarized in Table HU We ob- 
serve a general systematic shift more or less significantly 
below unity. These reevaluations unveil a new reactor 
antineutrino anomaly, clearly illustrated in Figure [SJ but 
still to be explained. In order to quantify the statistical 
significance of the anomaly we can compute the weighted 
average of the ratios of expected over predicted rates, for 
all short baseline reactor neutrino experiments (including 
their possible correlations). 

We consider the following experimental rate informa- 
tion: Bugey-4 and Rovno91, the three Bugey-3 experi- 
ments, the three Goesgen experiments and the ILL ex- 
periment, the three Krasnoyarsk experiments, the two 
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# 


result 


Det. type 


Tn (s) 


B D 


U3»p u 




a4i Pu 


old 


new 


err(%) 


corr(%) 


L(m) 


1 


Bueev-4 


J He+H 2 


888.7 


0.538 


0.328 


0.078 


0.056 


0.987 


0.942 


3.0 


3.0 


15 


2 


R0VN091 


3 He+H 2 


888.6 


0.614 


0.274 


0.074 


0.038 


0.985 


0.940 


3.9 


3.0 


18 


3 


Ruffev-3-I 


"Li-LS 


889 


0.538 


0.328 


0.078 


0.056 


0.988 


0.946 


4.8 


4.8 


15 


4 


Bugey-3-II 


6 Li-LS 


889 


0.538 


0.328 


0.078 


0.056 


0.994 


0.952 


4.9 


4.8 


40 


5 


Bugey-3-III 


6 Li-LS 


889 


0.538 


0.328 


0.078 


0.056 


0.915 


0.876 


14.1 


4.8 


95 


6 


Goesgen-I 


3 He+LS 


897 


0.620 


0.274 


0.074 


0.042 


1.018 


0.966 


6.5 


6.0 


38 


7 


Gocsgen-II 


3 He+LS 


897 


0.584 


0.298 


0.068 


0.050 


1.045 


0.992 


6.5 


6.0 


45 


8 


Goesgen-I I 


3 He+LS 


897 


0.543 


0.329 


0.070 


0.058 


0.975 


0.925 


7.6 


6.0 


65 


9 


ILL 


3 He+LS 


889 


~ 1 








0.832 


0.802 


9.5 


6.0 


9 


10 


Kr£Lsn. ] 


3 He+PE 


899 


~ 1 








1.013 


0.936 


5.8 


4.9 


33 


11 


Krasn. I] 


3 He+PE 


899 


~ 1 








1.031 


0.953 


20.3 


4.9 


92 


12 


\\ VflQTl 1 1 i 
1V1 doll. _L_L_L 


3 He+PE 


899 


~ 1 








0.989 


0.947 


4.9 


4.9 


57 




SRP T 

O I \ I 1 


Gd-LS 


OO ( 


~ 1 








D QR7 




^ 7 


3 7 


1 s 


14 


SRP II 


Gd-LS 


887 


~ 1 


— 


— 


— 


1.055 


1.018 


3.8 


3.7 


24 


15 


ROVN088-1I 


3 He+PE 


898.8 


0.607 


0.277 


0.074 


0.042 


0.969 


0.917 


6.9 


6.9 


18 


16 


ROVN088-2I 


3 He+PE 


898.8 


0.603 


0.276 


0.076 


0.045 


1.001 


0.948 


6.9 


6.9 


18 


17 


ROVN088-1S 


Gd-LS 


898.8 


0.606 


0.277 


0.074 


0.043 


1.026 


0.972 


7.8 


7.2 


18 


IS 


ROVN088-2S 


Gd-LS 


898.8 


0.557 


0.313 


0.076 


0.054 


1.013 


0.959 


7.8 


7.2 


25 


19 


ROVN088-3S 


Gd-LS 


898.8 


0.606 


0.274 


0.074 


0.046 


0.990 


0.938 


7.2 


7.2 


18 



TABLE II. iVobs/iVpred ratios based on old and new spectra. Off-equilibrium corrections have been applied when justified. The 
err column is the total error published by the collaborations including the error on Stot, the corr column is the part of the error 
correlated among multiple-baseline experiments, or experiments using the same detector. This table is used to construct the 
covariance matrix used in Eq. 1101 



Savannah River results (SRP), and the five Rovno88 ex- 
periments. R is the corresponding vector of 19 ratios of 
observed to predicted event rates. We assume a 2.0% sys- 
tematic uncertainty fully correlated among all 19 ratios. 
This choice is motivated by the common normalization 
uncertainty of the corresponding beta-spectra measured 
in (2l| - [23 |. We considered the ratios and relative errors 
gathered in Table [TTJ In order to account for the po- 
tential experimental correlations, we fully correlated the 
experimental errors of Bugey-4 and Rovno91, of the three 
Goesgen and the ILL experiments, the three Krasnoyarsk 
experiments, the five Rovno88 experiments, and the two 
SRP results. We also fully correlated the Rovno88 
(II and 21) results with Rovno91, and we added an arbi- 
trary 50% correlation between the Rovno88 (II and 21) 
and the Bugey-4 measurement. We motivated these lat- 
est correlations by the use of similar or identical integral 
detectors. We stress here that in this publication we 
only used the error budget published by the collabora- 
tions, without any change. We then obtain the covari- 
ance matrix W of the ratios. In Figure Q] we show the 
corresponding correlation matrix with labels detailed in 
Table HI 

All the above-mentioned experiments published ratios 
of measured to expected event rates. While the rates 
themselves can be considered to follow Gaussian distri- 
butions, the ratio of two Gaussian variables does not. 
This may lead to overestimating the statistical signifi- 
cance of the deviation of the ratios from unity, as was 
pointed out by James in response to a paper by Reines 
et al. 66]. In order to check our results, we performed 




FIG. 1. Correlation matrix of 19 measurements of reactor 
neutrino experiments operating at short baselines. Experi- 
ment labels are given in Table ITT1 



simple Monte-Carlo simulations: each experiment's ob- 
served and expected rates were simulated within their 
quoted errors, including the correlations between exper- 
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iments shown in the covariance matrix W. 

We then calculated the best-fit average ratio r bcst 
for each simple Monte-Carlo experiment by minimizing 
the x 2 function with respect to r: 

(r- lt) 7 V~ 1 (r- it) , (10) 

and obtained the distribution of rb os t (we note \i its 
mean in the following paragraph). We found that 
it is almost Gaussian, but with slightly longer tails, 
which we decided to take into account in our calcula- 
tions (in contours that appear later in this article we 
enlarged the error bars). If we only consider experi- 
ments with baselines < 100 m to get rid of a possible 
(6*i3, Amf^) driven oscillation effect at Palo Verde or 
CHOOZ, with the old antineutrino spectra the mean ra- 
tio is ^=0.976±0.024, and the fraction of simple Monte- 
Carlo experiments with r > 1 is 17.1% (— 0.95cr from 
expectation). With the new antineutrino spectra, we ob- 
tain ^=0.943±0.023, and the fraction of simple Monte- 
Carlo experiments with r > 1 is 1.3%, corresponding to 
a —2.2(7 effect (while a simple calculation assuming nor- 
mality would lead to — 2.4er). Clearly the new spectra 
induce a statistically significant deviation from the ex- 
pectation. In the following we define an experimental 
cross section af i0 = 0.943 x cr P rcd ' now 10 _ 43 cm 2/ fission . 
With the new antineutrino spectra, we observe that for 
the data sample the minimum x 2 is Xmin data = 19.6. 
The fraction of simple Monte-Carlo experiments with 
Xmin < X m in, data ^ s 25%, showing that the distribution 

of experimental ratios in R around the mean value is 
representative given the correlations. 

We will now discuss the possible explanations of this 
deviation from unity: an erroneous prediction of the an- 
tineutrino flux from the reactors, or a correlated artifact 
in the experiments, or a real physical effect if both pre- 
vious cases are excluded. 

Due to the importance of the antineutrino rate increase 
we suggest that independent nuclear physics groups 
should perform similar computations. We also consider 
that new measurements of the electron spectra of irra- 
diated fissile isotopes would help clarifying the anomaly. 
All cross sections of reactions used for the absolute and 
relative normalizations of the ILL electron spectra have 
been checked and found in agreement with the published 
values within error bars. A more complete discussion on 
the evaluation of the normalization of reactor antineu- 
trino spectra based on the ab initio method will be pub- 
lished later in [5l| . 

Assuming the correctness of CT P red ' new ^he anomaly 
could be explained by a common bias in all reactor neu- 
trino experiments. The measurements used one of two 
techniques, scintillator counters and integral detectors. 
The Bugey-3 experiment [|| used 3 identical detection 
modules, each of 600 liters, filled with 6 Li-loaded liquid 
scintillator. Bugey-3 recorded 120,000 neutrino interac- 
tions. The Bugey-4 experiment [|| used the Rovno91 Q 
integral type detector, but increasing the antineutrino 



rate by a factor of three. A similar detector was used 
for two integral measurements Rovno88 II and 21 Q- 
In such detectors, based on 3 He-filled counters, the to- 
tal number of antineutrino interactions is measured by 
detecting only the neutrons from reaction Eq. [3] The 
Goesgen experiment [B| used a detector nearly identical 
to the one used for the ILL neutrino experiment 0, but 
with the additional feature of position sensitivity. More 
than 10,000 neutrino events were recorded at the three 
detector locations. The detector contained liquid scintil- 
lator surrounded by 3 He-filled wire chambers for neutron 
detection. Both the positron and the neutron from re- 
action Eq. [3] were detected. Krasnoyarsk [1, Q used an 
integral detector with a scintillation section. The Savan- 
nah River experiments considered in this article used a 
scintillator counter |ll| . 

Neutrons were tagged either by their capture in metal- 
loaded scintillator, or in proportional counters, thus lead- 
ing to two distinct systematics. As far as the neutron de- 
tection efficiency calibration is concerned, we note that 
different types of radioactive sources emitting MeV or 
sub-MeV neutrons were used (Am-Be, 252 Cf, Sb-Pu, Pu- 
Be). 

It should be mentioned that the Krasnoyarsk, ILL, 
and SRP experiments operated with nuclear fuel such 
that the difference between the real antineutrino spec- 
trum and that of pure 235 U was less than 1.5%. They 
reported similar deficits to those observed at other reac- 
tors operating with a mixed fuel of 235 U, 239 Pu, 241 Pu 
and 238 U. Hence the anomaly cannot be associated with 
a single fissile isotope. 

All the elements discussed above argue against a trivial 
bias in the experiments, but a detailed analysis of the 
most sensitive of them, involving experts, would certainly 
improve the quantification of the anomaly. 

As discussed in Section UH in a near future the reactor 
rate anomaly significance might evolve due to the reevalu- 
ation of the cross section prefactor k to 0.3961 10 -43 cm 2 . 
The averaged ratio would shift to 0.938±0.023, leading 
to a deviation from unity at 99.2% C.L. 

The third kind of possible explanation of the anomaly 
based on a real physical effect will be detailed in Sec- 
tion El 



B. Shape information 

In this Section we re-analyze the Bugey-3 and ILL 
shape information, based on the published data @, H[. 
We will use this information for our combined analysis 
described in the next section. 



1. Bugey-3 

Based on the analysis of the shape of their energy spec- 
tra at different source-detector distances [1, Q , the Goes- 
gen and Bugey-3 measurements exclude oscillations such 
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that 0.06 < Am 2 < 1 eV 2 for sin 2 (2(9) > 0.05. 

For further analysis we used Bugey-3's 40 m/15 m ratio 
data from [3| as it provides the best limit. We followed 
the steps outlined in Q, building the following \ 2 func- 
tion: 

where R l ohs are our simulated data for Bugey from our 
tuned simulation and R[ h are our Monte-Carlo expecta- 
tion in each bin. The Oi are the errors reported by the 
Bugey-3 collaboration, and a is a systematic parameter 
accounting for the o~ a = 2% uncertainty on the relative 
normalization at 40 m and 15 m. Figure 16 of [4] shows 
the 90% C.L. exclusion contour from a raster scan anal- 
ysis with this estimator. As can be seen in Figure [2] we 
adequately reproduce Bugey-3's results. 
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We reanalyzed the data of Ref. [36| in order to in- 
clude the ILL shape-only information in our analysis of 
the reactor antineutrino anomaly. We built the following 
shape-only estimator 

2 _ 'sT' 6 f ( 1 + a ) R th ~ R obs \ 2 / 10 n 

AILL, shape — 2-^1 \ / ^ > 

i=l ^ 1 ' 

where i?o bs are the measured ratios in each energy bin, 
and R l th are our Monte-Carlo expectation in each bin. a 
is a free parameter in the fit, which renders this estimator 
completely insensitive to any normalization information. 
It is therefore only sensitive to the shape of the distribu- 
tion. 

Gi is the total error in each bin: we added in quadra- 
ture the statistical error and a systematic error of 11%. It 
was difficult to extract the magnitude of this shape-only 
systematic error from published information. We com- 
bined our Xill shape estimator with the rate-only estima- 
tor, and verified that with this value of the systematic, we 
could reproduce the ILL contours of both papers: when 
using i?iLL.80 we reproduced the contour in Figure 14 of 
Ref. J2|, and when using -Rill, 95 we reproduced that of 
Ref. |36| which excludes the no-oscillation hypothesis at 
2cr. A systematic error of 11% is consistent with Fig- 
ure 13 of Ref. [2]. It is also a conservative treatment of 
the shape-only information: with such an error the data 
are compatible with the no-oscillation hypothesis at la. 
Figure [3] shows the data points from 36 , along with our 
best fit from the shape-only estimator (|Am 2 | ~ 2.3 eV 2 , 
sin 2 (2(9) ~ 0.24), and the no-oscillation line for compari- 
son. 



FIG. 2. 90% C.L. exclusion domains obtained in the Am 2 - 
sin 2 (2#) plane from a raster scan of Bugey-3's data. Our 
result (continuous line) is in good agreement with the original 
result from 0] (dashed line), excluding oscillations such that 
0.06 < Am 2 < 1 eV 2 for sin 2 (20) > 0.05. 



We already note here that in Section |Vl when combin- 
ing Bugey-3's shape information with other results, we 
use the estimator from equation 111! but we perform a 
standard global scan (i.e. the minimization is performed 
over the entire Am 2 , sin 2 (29) plane). We only perform a 
raster scan here to show our agreement with the collab- 
oration's original analysis. 




3 4 

Energy (MeV) 



2. ILL 

As already noted in Ref. [36[ the data from ILL showed 
a spectral deformation compatible with an oscillation 
pattern in their measured over predicted events ratio. 
It should be mentioned that the parameters best fit- 
ting the data reported by the authors of Ref. [36| were 
Am 2 = 2.2 eV 2 and sin 2 (2(9) = 0.3. 



FIG. 3. Ratio of measured to expected positron energy spec- 
tra of the ILL neutrino experiment (data points extracted 
from [3^]). We show the best fit line with oscillations, along 
with the no-oscillation line, from our shape-only fit. The short 
error bars are statistical and the longer ones include the 11% 
systematic error described in the text. 
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V. THE FOURTH NEUTRINO HYPOTHESIS 

In this section we discuss the compatibility of the reac- 
tor antineutrino anomaly with the existence of a fourth 
non-standard neutrino, corresponding in the flavor basis 
to the existence of a sterile neutrino v s (see [l[ and refer- 
ences therein). The motivation is the explanation of the 
antineutrino deficit by an oscillation of electron neutrinos 
into a new neutrino state with a large Ara 2 ew value. 

For simplicity we restrict our analysis to the 3+1 four- 
neutrino scheme in which there is a group of three ac- 
tive neutrino masses separated from an isolated neutrino 
mass, such that |A?7i 2 ow | 3> 10 -2 eV 2 . The latter would 
be responsible for very short baseline reactor neutrino 
oscillations. For energies above the inverse beta decay 
threshold and baselines below 2 km, we adopt the ap- 
proximated oscillation formula from Ref . 52] : 



P, 



1 - cos 4 ncw sin 2 (26»i 3 ) sin 2 ' Am ^ L 



*E . 



sin 2 (20 now ) sin 2 



Am 



4Ek 



(13) 



1 dof Ax profile 



90.00 % 
95.00 % 
99.00 % 
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FIG. 4. Allowed regions in the sin 2 (2# ncw ) — Am^ cv plane 
obtained from the fit of the reactor neutrino data, without 
ILL-shape information, but with the stringent oscillation con- 
straint of Bugey-3 based on the 40 m/15 m ratios to the 3+1 
neutrino hypothesis, with sin 2 (26*13) = 0. The best-fit point 
is indicated by a star. 



In our analyses the well known solar driven oscillation ef- 
fects are negligible. The contribution of the atmospheric 
driven oscillation is negligible at distances < 100 m. It 
is worth noting that we are not sensitive to any sterile 
neutrino mass hierarchy in the mass range considered. 

The ILL experiment may have seen a hint of oscil- 
lation in their measured positron energy spectrum (see 
Section IIIIB|) . but Bugey-3 's results do not point to 
any significant spectral distortion more than 15 m away 
from the antineutrino source. Hence, in a first approx- 
imation, hypothetical oscillations could be seen as an 
energy-independent suppression of the u e rate by a fac- 
tor of ~ sin 2 (26 ) nC w,R), thus leading to Am 2 ew R > 1 eV 2 
and accounting for Bugey-3 and Goesgen shape analy- 
ses [3, H[ . Considering the weighted averaged of all reac- 
tor experiments we get an estimate of the mixing angle, 
sin 2 (20 now , R ) -0.115. 

Let us now fit the sterile neutrino hypothesis to the 
data (baselines below 2 km) by minimizing the least- 
squares function 



(P ee -lt) ^(Pea-it) 



(14) 



assuming sin 2 (26 ) i 3 ) = 0. 

Figure @] provides the results of the fit in the 
sin 2 (26 J ncw ) — Am 2 cw plane, including the reactor ex- 
periment rate information as well as the 40 m/15 m 
Bugey-3 spectral shape constraint presented in Figure [2] 
The latter leads to stringent oscillation constraints for 
|Ato 2 cwR | < 1 eV 2 , since no spectral distortion was 
observed in the ratio 40 m/15 m. The fit to the 
data indicates that |Am 2 cwR | > 0.2 eV 2 (95%) and 

sin (26* ne w,R,) — 0.1. Note that if we include a non-zero 
value of sm 2 (26»i 3 ), as large as the 90% C.L. CHOOZ 



bound 13], the contours presented in Figure [4] are only 
marginally affected. 

Ignoring the reactor anomaly leads to an ambiguous in- 
terpretation of the —1-2 km baseline experiment results 
constraining #13. As an example we compare the follow- 
ing two hypotheses: the 3 neutrino mixing scenario with 
sin 2 (26'i3) = 0.06 for instance and the 3+1 neutrino mix- 
ing scenario, taking sin 2 (26*13) = 0, sin 2 (20 nC w,R.) = 0.115 
and I Am 2 ew R | > 1 eV 2 (best values of the combined 
fit shown on Figure E|) . Figure [5] displays the data to- 
gether with the two hypotheses. The red line repre- 
sents the three active neutrino hypothesis, considering 
reactor neutrino rates until the CHOOZ baseline. The 
p-value of the 3 neutrino mixing hypothesis described 
above is 17% (26.9/21 degrees of freedom). The blue 
line illustrates the 3+1 neutrino mixing scenario with a 
large |Am 2 ew R | = 2.4 eV 2 , leading to a p-value of 55% 
(19.4/21 degrees of freedom). We note a tension between 
the data and the three active neutrino flavor hypothesis, 
again illustrating the anomaly. Clear preference is given 
to the 3+1 neutrino hypothesis since A\ 2 > 7. A simi- 
lar computation with the former reference spectra would 
lead to a p-value of 34% in both cases. This illustrates the 
impact of our reanalysis with the new antineutrino spec- 
tra, updated neutron mean lifetime, and off-equilibrium 
effects. 

At this stage it is tempting to consider the previously 
noted anomalies affecting other short baseline electron 
neutrino experiments Gallex, Sage and MiniBooNE, re- 
viewed in Ref. [56| . Our goal is to quantify the compati- 
bility with those anomalies. 

We first reanalyzed the Gallex and Sage calibra- 
tion runs with 51 Cr and 37 Ar radioactive sources emit- 
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FIG. 5. Illustration of the short baseline reactor antineutrino anomaly. The experimental results are compared to the prediction 
without oscillation, taking into account the new antineutrino spectra, the corrections of the neutron mean lifetime, and the 
off-equilibrium effects. Published experimental errors and antineutrino spectra errors are added in quadrature. The mean 
averaged ratio including possible correlations is 0.943 ± 0.023. The red line shows a possible 3 active neutrino mixing solution, 
with sin 2 (2#i3) = 0.06. The blue line displays a solution including a new neutrino mass state, such as |Am 2 ew R | 3> 1 eV 2 and 
sin 2 (2# ncWi R,) = 0.12 (for illustration purpose only). 



ting ~ 1 MeV electron neutrinos. l57f, following the 
methodology developed in Ref. [H, 158]. However we 
decided to include possible correlations between these 
four measurements in this present work. Details are 
given in Appendix [Bj This has the effect of being 
slightly more conservative, with the no-oscillation hy- 
pothesis disfavored at 97.7% C.L., instead of 98% C.L. 
in Ref. [56| . Gallex and Sage observed an average deficit 
of Rq = 0.86 ± 0.06 (la). Considering the hypothesis of 
v e disappearance caused by short baseline oscillations we 
used Eq. (flU]) . neglecting the Am^ driven oscillations 
because of the very short baselines of order 1 meter. Fit- 
ting the data leads to |Am£ eW]Q | > 0.3 eV 2 (95%) and 

sin 2 (26' n0Wi G) ~ 0.26. Combining the reactor antineu- 
trino anomaly with the gallium anomaly gives a good fit 
to the data and disfavors the no-oscillation hypothesis at 
99.7% C.L. Allowed regions in the sin 2 (26 , ne w) — Am 2 cw 
plane are displayed in Figure [5] (left). The associated 
best-fit parameters are |Am 2 cw R&G | > 1.5 eV 2 (95%) 
and sin 2 (26> ncWiR&G ) ~ 0.12. 

We then reanalyzed the MiniBooNE electron neutrino 
excess assuming the very short baseline neutrino os- 
cillation explanation of Ref. (56|. Details of our re- 
production of the latter analysis are provided in Ap- 
pendix[B] The best fit values are |Am 2 ew MB | = 1.9 eV 2 

and sin 2 (2# ncW: MB) ~ 0.2, but are not significant at 
95% C.L. The no-oscillation hypothesis is only disfa- 
vored at the level of 72.4% C.L., less significant than 
the reactor and gallium anomalies. Combining the re- 
actor antineutrino anomaly with our MiniBooNE re- 



Experiment (s) 


sin 2 (2# ncw ) 


|Am;U| (eV) 


C.L. (%) 


Reactors (no ILL-S,R*) 


0.02-0.20 


> 0.40 


96.5 


Gallium (G) 


> 0.06 


> 0.13 


96.1 


MiniBooNE (M) 






72.4 


ILL-S 






68.1 


R* + G 


0.05-0.22 


> 1.45 


99.7 


R* + M 


0.04-0.20 


> 1.45 


97.6 


R* + ILL-S 


0.02-0.21 


> 0.23 


95.3 


All 


0.06-0.22 


> 1.5 


99.8 



TABLE III. Best fit parameter intervals or limits at 95% C.L. 
for sin 2 (2# ncw ) and |Am 2 cw | parameters, and significance of 
the sterile neutrino oscillation hypothesis in %, for different 
combinations of the reactor experimental rates only (R*), the 
ILL-energy spectrum information (ILL-S), the gallium experi- 
ments (G), and MiniBooNE-^ (M) re-analysis of Ref. [5(|. We 
quantify the difference between the sin 2 (2# ncw ) constraints 
obtained from the reactor and gallium results. Following pre- 
scription of Ref. [77|, the parameter goodness-of-fit is 27.0%, 
indicating reasonable agreement between the neutrino and an- 
tineutrino data sets (see Appendix |B)) . 



analysis leads to a good fit with the sterile neutrino 
hypothesis and disfavors the absence of oscillations at 
98.5% C.L., dominated by the reactor experiments data. 
Allowed regions in the sin 2 (2# ncw ) — Am 2 ow plane are 
displayed in Figure [S] (right). The associated best- 
fit parameters are |Am 2 ewR&MB | > 0.4 eV 2 (95%) and 

sin 2 (20 new 3&MB) ~ 0.1. 
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1 dof Ax 2 profile 1 dof Ax 2 profile 




FIG. 6. Allowed regions in the sin 2 (2t5 ncw ) — Am„ ew plane obtained from the fit of the reactor neutrino data to the 3+1 neutrino 
hypothesis, with sin 2 (2#i3) = 0. The left panel is the combination of the reactors and the gallium experiment calibration results 
with 51 Cr and 37 Ar radioactive sources. The right panel is the combination of the reactors and our reanalysis of the MiniBooNE 
data following the method of Ref. 56|. In both cases the ILL energy spectrum information is not included. 



Our ILL re-analysis, including only the en- 
ergy spectrum shape, leads to the allowed regions 
in the sin 2 (2(9 new ) — Ato 2 gw plane presented in 
Figure [7] We notice a hint of neutrino oscil- 
lations SUCh that |ATOncw,ILL-shapel > 1 ^ and 

sin 2 (26> noW! T LL _ shapc ) ~ 0.2, in agreement with our 
fourth neutrino hypothesis, but still compatible with the 
absence of oscillations at the la level. Figure [3] is our 
reproduction of the illustration 3 of Ref. Q ; we superim- 
posed the oscillation pattern that would be induced by 
neutrino oscillations at our best fit (combined analysis). 
The ILL positron spectrum is thus in agreement with 
the oscillation parameters found independently in our 
re-analyses, mainly based on rate information. Because 
of the differences in the systematic effects in the rate 
and shape analyses, this coincidence is in favor of a true 
physical effect rather than an experimental anomaly. As 
a cross check we performed a Monte- Carlo simulation of 
the ILL and Bugey-3 experiments, including the finite 
spatial extension of the nuclear reactors and the ILL and 
Bugey-3 detectors. We found that the small dimensions 
of the ILL nuclear core lead to small corrections of the 
oscillation pattern imprinted on the positron spectrum. 
However the large extension of the Bugey nuclear core is 
sufficient to wash out most of the oscillation pattern at 
15 m. This explains the absence of shape distortion in 
the Bugey-3 experiment. 

Table IIIII summarizes all the results of our fits of reac- 
tor, gallium, and MiniBooNE- v data to the sterile neu- 
trino oscillation hypothesis. We observe that all the data 
sets taken separately are very consistent with one an- 
other, pointing to very similar oscillation parameters. We 
thus performed a global fit to all available data. 

The no-oscillation hypothesis is disfavored at 
99.8% C.L. The significance is dominated by the 
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FIG. 7. Allowed regions in the sin 2 (2t9 ncw ) — Am 2 ew plane 
obtained from a fit of the ILL energy spectrum shap e only. 
The best fit value reported by the authors of Ref. [3(| is very 
close to our best fit, at |Am 2 ew | ~2 eV 2 , but it is worth noting 
its poor statistical significance, compatible with the absence 
of oscillations at the la level. The best-fit point is indicated 
by a star. 



gallium and reactor data. Allowed regions in the 
sin 2 (2(9 ncw ) — Am 2 0W plane are displayed in Figure |Sl 
together with the marginal A^ 2 profiles for |Ato 2 0W | and 
sin 2 (2(9 ncw ). The combined fit leads to the following con- 
straints on oscillation parameters: |Am 2 cw | > 1.5 eV 2 
(95% C.L.) and sin 2 (2(9 ncw ) = 0.14 ± 0.08 (95% C.L.). 
An embryo of possible consequences of this result will 
be discussed in Section [Villi 
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FIG. 8. Allowed regions in the sin 2 (2# ncw ) — AmJ„ plane from the combination of reactor neutrino experiments, Gallex and 
Sage calibration sources experiments, MiniBooNE reanalysis of Ref. [5f| , and the ILL-energy spectrum distortion. The data are 
well fitted by the 3+1 neutrino hypothesis, while the no-oscillation hypothesis is disfavored at 99.8% C.L. The marginal A\ 2 
profiles for |Am 2 cw | and sin 2 (26» ncw ) (1 dof) lead to the constraints, !Am 2 ew | > 1.5 eV 2 (95% C.L.) and sin 2 (26> ne w) = 0.14±0.08 
(95% C.L.). 



VI. REVISION OF THE CONSTRAINTS ON 13 

In this section we discuss the impact of our revision 
of the results of short baseline (< 100 m) reactor exper- 
iments on the constraints on the #13 mixing angle from 
CHOOZ [H, KamLAND [l7j], and their combination. 



The results depend on the choice of the cross section per 
fission used to normalize the predicted antineutrino rate, 
as well as on the neutrino oscillation scheme used in the 
analysis, i.e. whether to include a non-standard neutrino 
or not. Assuming the correctness of the neutrino exper- 
iments (detected rate), we show that constraints on #13 
could be derived by using cr^ re ,ano and the three-active 
neutrino scenario. 



A. CHOOZ 

We first reproduced CHOOZ's background-subtracted 
results and Monte-Carlo using Ref. [13] • We used the 
converted neutrino spectra from ILL electron data [22I 
|40L l4lj . using the parameterization from Ref. 47 1, and 
the V-A cross section from [29|, HH ^normalized to the 
Bugey-4 value as was done in Ref. Long-lived fis- 

sion product contributions were also taken into account 
(at 300 days of irradiation). As in Ref. [l^ . we accounted 



for a a a = 2.7% systematic error on the overall normal- 
ization and a 1.1% error on the energy scale, along with 
a bin-to-bin uncorrelated uncertainty describing the un- 
certainty of the neutrino spectrum. Using the original 
ILL neutrino converted spectra, with the a 
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section we could roughly reproduce the exclusion limit 
from Ref. [l3[ (black dashed line on Figure to ac- 
ceptable precision (gray line on Figure in the relevant 
Am| x range, although we could not fully imitate every 
step in the published analysis. In particular we did not 
use the same statistical treatment, using only a global 
scan while the CHOOZ collaboration used the Feldman- 
Cousins unified approach to extract the confidence inter- 
val. 

Let us now revisit the CHOOZ results using the new 
antineutrino spectra in the three active neutrino frame- 
work. We perform the same analysis as before, but 
computing the expected number of events based on 
ore , new an( j ^ g error CTq was ^us increased to 3.3% 

to account for this effect. This leads to a new exclu- 
sion limit, sin 2 (26»i3) < 0.22 (1 dof) at 90% C.L. for 
Awijj = 2.4 10 -3 eV 2 . The corresponding exclusion con- 
tour is displayed as the blue line on Figure [SI In this ap- 
proach we make use of a P ied < new an d attribute the slight 
deficit of v e at CHOOZ to a 6*13 driven oscillation, such 



that sin^(26»i 3 ) = 0.11 



+0.07 
-0.08- 



However this result only 
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holds within the three active neutrino mixing framework, 
under the hypothesis of no-oscillation at very short base- 
lines. But in Section IIVI we pointed out an anomaly 
affecting the short baseline reactor neutrino experiment 
results. In the scenario discussed above with a new neu- 
trino (Am 2 ow , 6 ne w) we generate oscillations driven by a 
large |Am 2 ew | ^> Am^. In this case CHOOZ measured 
a combination of the oscillation driven by (Am 2 ew , # new ) 
and (Amjj, #13). A comprehensive 3+1 neutrino analysis 
would thus be mandatory to constrain # 13 using the new 
normalization based on a P lcd ' ncw ^ Du |; we w [\\ see nex t 
that another approach may be used. 
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FIG. 9. Our simulation (gray line) roughly reproduces the 
published CHOOZ exclusion limit [13(] represented by the 
black dashed line. We obtain the blue contour by normalizing 
the expected no-oscillation antineutrino rate to <jP rcd > ncw an d 
analyzing the data using the three active neutrino oscillations 
scheme. We then normalize the expected antineutrino event 
rate to cr| no to account for possible neutrino oscillations at 
very short baselines, leading to the green contour. This pro- 
vides our new constraint on sin 2 (2613), accounting for the 
reactor antineutrino anomaly. 

As a matter of fact the CHOOZ experiment used 
<7j Ugey instead of a P rod,now . The collaboration justified 

making use of ugey because of the good agreement with 

er? re '° . With our revision of <r^ red this justification does 
not hold anymore in the 3 neutrino oscillation framework. 
But it turns out that if we assume non-standard neu- 
trino oscillations in the large |Am 2 ow | ^ 0.1 eV 2 regime 
(no spectral distortion), the normalization of CHOOZ 
using o-? ugey ~ Oy no leads to an estimator of sin 2 (2#i 3 ) 

which eliminates a possible degeneracy with sin 2 (2# new ). 
CHOOZ's strategy was indeed to absorb possible errors 
in rate predictions, but this methodology holds as well for 
constraining #13 on top of an additional, short baseline, 
averaged oscillation. 

Thus a pragmatic approach for constraining #13 from 
CHOOZ data is the use of the weighted average of the 
measured cross sections of all experiments < 100 m, cr^ no . 



This leads to a new exclusion limit, sin 2 (2013) < 0.10 
(1 dof) at 90% C.L. for Am^ = 2.4 10" 3 eV 2 , slightly 
lower than CHOOZ's published value (T3 |. 

B. KamLAND 

In this section we reevaluate the KamLAND constraint 
on 0i3 in light of the reactor antineutrino anomaly. 

Unlike in the original CHOOZ analysis from 13], 
the KamLAND D e energy spectrum is calculated using 
nred.old j-, ase d on Ref. Ho|. But the spectral uncer- 
tainty is evaluated from Ref . Q , including off-equilibrium 
effects. In this case the propagation of the Bugey-4 error 
on <7? red only has a marginal impact since it is not among 
the dominant systematics. 

We first reproduced KamLAND 's results and Monte- 
Carlo to a very good accuracy exploiting Ref. (l7j as 
well as publicly available information [4a]. We used 
the parameterization of Ref. [47| for Stotj and included 
off-equilibrium corrections according to Ref. [20]. Spe- 
cial care was taken to include all published backgrounds, 
known antineutrino sources, especially Korean power re- 
actors, and geoneutrinos. We tuned our simulation to 
reproduce KamLAND's Monte-Carlo with and without 
neutrino oscillations (adjusting each effective power of 
the Japanese power plants). In both cases our simula- 
tion agrees with that of KamLAND to better than 1% in 
the 1.5-6.5 MeV range. 

Confidence levels in the (tan 2 #12, sin 2 #13) plane were 
obtained by minimizing the generic \ 2 function: 

^ = £(^f^)V»^ (15) 

where the Yi are our simulated data tuned to Kam- 
LAND's Monte-Carlo [13] ■ The free parameters are a 
($12, #13, Anijj and geo-is rate), and the nuisance param- 
eters, P, which are the systematics quoted in Ref. p7| . 
The Ni(a, (3) are our simulation model for all the free pa- 
rameters. The W~ 1 matrix contains all the systematic 
uncertainties quoted in Ref. \47\ . 

Figure [TO] demonstrates that we could reproduce the 
best fit and the confidence contours of KamLAND [l7| 
with good accuracy. With the original reactor neu- 
trino spectra Stotj we obtain a good agreement with 
KamLAND's published results [17] . Changing the ref- 
erence spectra S to t according to Ref. [2(|, the best fit 
values and uncertainties on tan 2 #12 and Am^ are un- 
affected in a three active neutrino, oscillations frame- 
work, but the sin 2 (2^13) central value is shifted upwards 
to sin 2 (26'i3) = 0.2+0^3, consistent with zero at 1.5a-. 
The new 90% C.L. limit on sin 2 (26*13) would therefore 
increase to 0.41 (1 dof), marginalizing over the other fit 
parameters. 

However, as in the case of CHOOZ discussed above, 
these results do not take into account the reactor an- 
tineutrino anomaly at short baselines. In the 3+1 neu- 
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FIG. 10. Our simulation (gray line, gray disk) reproduces the 
published 95% C.L. KamLAND inclusion contour [TtJ repre- 
sented by the black dashed line (black circle). We obtain the 
blue contour (star) by normalizing the expected no-oscillation 
antineutrino rate to cr P red > new a nd analyzing the data using the 
three active neutrino oscillations scheme. We then normalize 
the expected antineutrino event rate using af 10 to account 
for possible neutrino oscillations at very short baselines, lead- 
ing to the green contour (diamond). This provides our new 
constraints in the tan 2 #i2-sin 2 613 plane, accounting for the 
reactor antineutrino anomaly. We note that the slight tension 
between the sin 2 813 best-fit values of solar (red square) and 
reactor data is reduced in our scenario. 



trino oscillation framework, the whole effect induced by 
the normalization shift is absorbed by the new oscillation 
at very short baselines, driven by Am 2 0W . Normalizing 
KamLAND 's data to aJ no leads to an estimator revising 
the constraint on #13 using the three-active neutrino os- 
cillation framework, such that sin 2 (26' 1 3) < 0.31 (1 dof) 
at 90% C.L. for Am 2 31 = 2.4 10~ 3 eV 2 . 



C. Combining CHOOZ and KamLAND 

In sections IVI Al and IVIBI we separately revisited the 
CHOOZ and KamLAND results. Here, we first volun- 
tarily ignore the reactor antineutrino anomaly at short 
baselines. We thus normalize both CHOOZ and Kam- 
LAND with a P Yed ' new . Combining these two results leads 

to a best-fit value of sin 2 (26>i 3 ) = 0.13lo;Qg, barely consis- 
tent with a null value of 6*13. The left panel of Figure [TT1 
shows the Av 2 profiles for CHOOZ (green), KamLAND 
(blue), and for their combination (red), as a function 
of sin 2 (2f?i3), marginalizing over the other parameters. 
In this scenario the revised limit is sin 2 (26'i3) < 0.23 
at 90% C.L. (1 dof), for the estimation of 6* 13 only, for 
Arnlj = 2.4 10" 3 eV 2 . These results arc consistent with 
analyses of various data sets indicating a 'hint' for a non- 
zero 0i3 [HI, [ill] . This is explained by the increase of the 
predicted event rate at CHOOZ and KamLAND, which 
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TABLE IV. Best-fit values for sin 2 (2^13) and la errors, con- 
sidering two possible normalizations with (J P rod ' new or <jj no . 
Limits are given for a 1 dof parameter estimation. 



is attributed to #13 only, in the three neutrino oscillations 
framework. 

But the previous estimation does not take into ac- 
count the reactor antineutrino anomaly. Let us now 
consider that both CHOOZ and KamLAND are nor- 
malized with Oj no , in order to constrain #13 using a 
three active neutrino oscillation framework. The right 
panel of Figure Ql] shows the A% 2 profiles for CHOOZ 
(green), KamLAND (blue), and for their combination 
(red), computed as before. In this scenario the best fit 
is at sin 2 (26li3) = 0.0±£o7- This lower best fit value is 
due to the fact that we predict fewer antineutrinos at 
CHOOZ and KamLAND, some of them having already 
oscillated into non-standard neutrinos. This leads to our 
revised constraint, sin 2 (26»i 3 ) < 0.11 at 90% C.L. (1 dof) 
for Amjj = 2.4 10 -3 eV 2 . Results are summarized in 
Table 13 

We conclude that the hint for a non zero value of 6*13, 
first presented in Ref. [18| . vanishes if one normalizes 
CHOOZ's and KamLAND's predicted antineutrino rates 
to crj n °. Furthermore we notice that the slight tension 
affecting the sin 2 6*13 best-fit values of solar and reac- 
tor data [l9| is reduced in our scenario involving a non- 
standard neutrino such that Am 2 ew 3> Am^. 

A full analysis of the data based on non-standard neu- 
trino oscillation models would lead to more accurate re- 
sults, beyond the scope of this article. 



D. Daya Bay, Double Chooz, and RENO 

Three new reactor neutrino oscillation experiments, 
Daya Bay [14], Double Chooz 0, and RENO [HI will 
soon be coming online to significantly improve our knowl- 
edge of reactor antineutrino rates at various baselines. 

Among these experiments Double Chooz is the only 
one to operate a first phase without a near detector. The 
evidence for a reactor antineutrino anomaly triggers the 
discussion of the choice of the cross section per fission to 
be used for constraining #13, cr P rod ' now or as we n as 
its associated uncertainty. 

Straightforwardly using <jP ied < new m the three neutrino 
oscillation framework to interpret a hypothetical antineu- 
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FIG. 11. Dependence of A\ 2 on sin 2 (2#i3) for CHOOZ (green), KamLAND (blue), and their combination (red). Comparing 
the data with the excess in expected event rate predicted using cr^ re ,new leads to an increase of the best-fit value of sin 2 (2#i3), 
in the three active neutrino oscillations framework (left). The normalization of the expected event rate with <jj n ° allows to 
absorb the reactor antineutrino anomaly observed at very short baseline, leading to our revised constraints on sin 2 (2^13). 



trino deficit at Double Chooz could lead to an overesti- 
mation of #13, which could fake a discovery. But if the 
sterile neutrino oscillation hypothesis discussed in this 
paper is proven, the use of -P rod ' now would be possible in 
active+sterile neutrino oscillation frameworks. In that 
case the 2.7% error budget of the antineutrino spectra 
should be used. Assuming that the average experimental 
results of reactor neutrino are correct, using <7j no has the 
advantage to absorb either a antineutrino flux miscalcu- 
lation, or a physical neutrino deficit at very short base- 
lines, leading to a conservative constraint on #13. The 
error budget could then be taken as the weighted stan- 
dard deviation of the short baseline experiments, 1.0%. 

This choice is however not relevant for experiments 
running with a multi-detector configuration, since they 
absorb part of the uncertainty of the reactor antineu- 
trino fluxes, depending on the setup. Even in the hy- 
pothetical case of antineutrino oscillation at very short 
baselines the sensitivities of Daya Bay, Double Chooz, 
and RENO should be marginally affected because of the 
large allowed values of Am^ ow . However this statement 
should be checked based on detailed simulations. 

In the non-standard neutrino hypothesis the discovery 
of a shape distortion in the positron energy spectrum of 
the far detectors may be determinant in disentangling #13 
from 6» now . 



VII. TESTING THE ANOMALY 

A. At reactors 

The presence of sterile neutrinos would leave its im- 
print on the signal at both the near and far detectors 
of forthcoming reactor neutrino experiments. Taking re- 



sults from multi-detectors and allowing for the possibility 
of sterile mixing angles, one can probe both #13 and the 
sterile mixing angle 9 new at Daya Bay, Double Chooz, 
and RENO [50] ■ Positron spectrum energy distortions 
should be deeply investigated at near detectors, as quoted 
in (53 |. In any case, measurements of the expected over 
predicted event rates at the near detectors will allow the 
probing of the reactor antineutrino anomaly, providing 
high statistics and high precision measurements at a few 
hundreds meters from the antineutrino sources. The an- 
tineutrino anomaly would be best tested by performing 
blind analyses for all near detector data. 

Further measurements at very short baselines below 
100 meters would be useful to confirm the MeV elec- 
tron antineutrino deficit. Currently no fundamental re- 
search program is underway to search for new oscilla- 
tion physics at reactors. However there is a worldwide 
program at short baselines for the purpose of nuclear 
non-proliferation, using antineutrinos as new IAEA safe- 
guards tools [Hj]. In this context the Nucifer experiment, 
located 7 meters away from the Osiris research reactor 
core in Saclay, will start its operation in 2011 [54j- Nu- 
cifer will thus have the possibility to test the anomaly. 
A rate-only analysis with a precision of a few percent 
may not be enough to provide a decisive improvement 
of the understanding of the anomaly. But a shape anal- 
ysis may provide enough information, depending of the 
energy resolution. 

We note here that the antineutrino non-proliferation 
program will not be affected by the antineutrino anomaly 
since relative antineutrino rates with respect to known 
thermal powers could be used to calibrate the experi- 
ments. 
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B. MegaCurie radioactive sources 

As mentioned above, radiochemical gallium experi- 
ments (Gallex and Sage) tested their experimental pro- 
cedure by exposing their gallium t arg et to MegaCurie 
neutrino sources using 51 Cr or 37 Ar [57|. The production 
and handling of such devices is thus well under control. 
Ref. |64| proposed to use a 51 Cr source inside two con- 
centric gallium tanks, whereas Ref. [6o| proposed to use 
a 51 Cr or 90 Sr source next to the Borexino detector. In 
liquid scintillating detectors 51 Cr or 37 Ar v^s are de- 
tected through neutrino-electron elastic scattering while 
90 Sr D e 's are detected through inverse beta decay. With 
the Am 2 ew values best fitting the sterile neutrino hypoth- 
esis, the deployment of a radioactive source at the center 
of an ultra-low background neutrino detector, such as 
Borexino, KamLAND, and SNO+, would allow both the 
testing of the v e deficit and the search for an oscillation 
pattern as function of the detector radius. These neu- 
trino sources emit quasi-mono-energetic neutrino lines of 
sub-MeV energies leading to a clear oscillation pattern 
at the range of a meter. In addition a 37 Ar source emits 
only 7-rays through second order processes and is there- 
fore easy to handle after its irradiation inside a breeder 
nuclear reactor. As an example, if a source of 1 MCi were 
inserted in the middle of a large detector like SNO+, it 
would provide a few hundred thousand interactions in 
the detector, of which a few ten thousands would deposit 
more than 500 keV. A < 15 cm accurate vertex recon- 
struction could allow to draw a simple and clear figure of 
the number of neutrino interactions as a function of ra- 
dius, directly testing the sterile neutrino oscillation pat- 
tern for Ato 2 < 5 eV 2 . 



VIII. CONCLUSIONS AND OUTLOOK 

The impact of the new spectra of reactor antineu- 
trino is extensively studied in this article. The increase 
of the expected antineutrino rate by about 3.5% com- 
bined with revised values of the antineutrino cross section 
significantly decreased the normalized ratio of observed 
to expected event rates in all previous reactor experi- 
ments performed over the last 30 years at distances below 
100 m Ellii. The average ratio is 0.943 ± 0.023, 
leading to the reactor antineutrino anomaly. This deficit 
could still be due to some unknown in the reactor physics, 
but we also analyze these revised results in terms of a 
suppression of the v e rate at short distance as could 
be expected from a sterile neutrino, beyond the stan- 



dard model, with a large |Am 2 GW | 3> 



lAmfjjJ. We note 
that hints of such results were already present at the 
ILL neutrino experiment in 1981 (36[. We also consid- 
ered that other neutrino experiments, MiniBooNE [59| 
and the gallium neutrino sources experiments [57|, ob- 
serve v e deficits at a similar L/E. These anomalies 
were comprehensively studied in Ref. [56| . It is impor- 
tant to note that these anomalies exist and are com- 



parable in both the neutrino and the antineutrino sec- 
tors. Furthermore it turns out that each experiment 
fitted separately leads to similar values of sin 2 (26> ncw ) 
and similar lower bounds for |Am 2 cw |, but without a 
strong significance. Hence, we combined in a global fit 
these results with our re-evaluation of the reactor ex- 
perimental results, taking into account the existing cor- 
relations. This leads to a solution for a new neutrino 
oscillation, such that |ATO 2 ew | > 1.5 eV 2 (95% C.L.) and 
sin 2 (26> now ) = 0.14±0.08 (95% C.L.), disfavoring the no- 
oscillation case at 99.8% C.L. This hypothesis should be 
checked against systematical effects, either in the predic- 
tion of the reactor antineutrino spectra or in the experi- 
mental results. 

We then revised the constraints on the #13 mixing an- 
gle obtained with the three active neutrino oscillation 
framework by CHOOZ [3 and KamLAND [l7|. We 
show that a D e deficit measured at any antineutrino ex- 
periment with a baseline greater than 1 km could be 
misinterpreted as a hint for a non-zero value of $13. Ac- 
counting for a possible | Am 2 cw |-driven sterile neutrino 
oscillation, we revise the constraint on the third mixing 
angle, such that sin 2 (26*13) < 0.11. Note that the Kam- 
LAND best fit of the solar neutrino parameters are un- 
changed. However, the combination of KamLAND and 
CHOOZ leads to a best fit value of sin 2 (26»i 3 ) = 0.0±g$, 
in good agreement with the best fit value extracted from 
the combined analysis of solar neutrino data. This re- 
laxes the tension between reactor and solar data recently 
reviewed in Ref. (l9j . 

If the existence of a fourth neutrino turns out to be 
verified we note that a 7% reduction of the total flux 
of solar neutrinos must be taken into account when 
confronting with the experimental results. In partic- 
ular the total neutrino v flux measured by SNO [68|, 
at 5.14±0.21 10 6 cm .s , is now in better agreement 
with the prediction of the (reduced) high-Z Sun Stan- 
dard Model model, at 5.25±0.9 10 6 cm^ls" 1 , and dis- 
favors the (reduced) low-Z one (4.09±0.9 10 6 cm _2 .s _1 ). 
We note that the high-Z model is also in good agree- 
ment with the data from helioseismology, contrarily to 
the low-Z model [6^ |. 

Assuming a hypothetical new neutrino i/ now heavier 
than the three active neutrinos we can briefly quantify 
its contribution to the effective masses searched for in 
/3-decay and neutrino-less /3/3-decay experiments, as per- 
formed in Ref. [56| . The Tritium /3-decay experiments 
have reported < 2 eV (95%) 67]. A fourth neutrino 
with a mass in the eV range should contribute more than 
the active neutrinos, which are expected to have sub- 
eV masses, to the signal searched for in direct detection 
beta decay experiments. We find that the contribution 
of the fourth neutrino state fitting the data analyzed in 
this article would contribute to mp for more than 0.2 eV 
(95% C.L.). This is within the sensitivity range of the 
forthcoming KATRIN experiment [7(3]. Assuming Majo- 
rana neutrinos, the contribution of the fourth state would 
be such that nipp > 0.02 eV (95% C.L.), which is above 
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the contribution of the three active neutrinos and thus 
disfavors possible cancellations with the three other ac- 
tive flavors [H], . This contribution is of interest for 
the forthcoming experiments (72[. 

Furthermore we would like to stress that the exis- 
tence of a fourth neutrino is slightly favored by some 
recent cosmological data. The effective number of neu- 
trino species fitted by WMAP and BAO observations [73[ 
is N off =4.34±0.87. An analysis of WMAP combined 
with the Atacama Cosmology Telescope data [zl] leads 
to Ncff =5.3±1.3, and an independent analysis of the 
WMAP 7 years data [iH provides N eff =4±l. Finally 
a recent analysis based on non-standard big bang nucle- 
osynthesis [76[ leads to N c ff=3.78±0.75. But the com- 
patibility of the sterile neutrino hypothesis exposed in 
this paper (|Am 2 ow | > 1.5 eV 2 , sin 2 (26> new ) = 0.14 ± 0.08 
at 95% C.L.) with cosmological models and data should 
be assessed, especially its contribution to the non- 
baryonic dark matter of the Universe, which may be non 
negligible. 

Finally, a clear experimental proof of the presence of 
this fourth non-standard neutrino becomes mandatory. 
This can be given by the imprint on the energy spectrum 
in a very short baseline reactor neutrino experiment or 
by a new neutrino source experiment in a detector with 
energy and spatial resolution. 
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Appendix A: Error propagation in reactor 
experiment analyses 

One of the main difficulties in these calculations is the 
propagation of the errors on the ILL data from [22|-[24j 
to the final result, i.e. the cross-section per fission for 
each of the isotopes. From our previous work in [20l | we 
propagated the error on the ILL data and the conversion 
procedure to the resulting neutrino spectra. However the 
binning and the energy range are then those of the ILL 
data (250 keV bins from 2 to 8 MeV for 235 U). 

These data are then fitted to an exponential- 
polynomial model, see for example [13]. This allows to 
use arbitrarily fine binning, which is necessary for the 
correct convolution with the inverse beta cross-section 
and also the oscillation probability. In (47| the authors 
outline a method to propagate the errors and correlations 
on the polynomial coefficients to the physical observable. 



However as they point out the correlation matrices ob- 
tained from the fit are rather unstable, most coefficients 
being strongly correlated or anti-correlated to each other. 
We found that this made error propagation very difficult. 

As an exercise, we performed the following simple 
Monte-Carlo simulation: using the original correlation 
matrix on the converted ILL spectra, we simulated a se- 
ries of ILL converted spectra, and fitted each of them us- 
ing the exponential-polynomial model from (47j . We then 
evaluated the fitted polynomials in each energy bin, yield- 
ing fitted Monte-Carlo neutrino spectra. We then com- 
puted the bin-to-bin correlation matrix of these Monte- 
Carlo spectra. 

An example of the result is shown in Figure I13[ along 
with the original correlation matrix in Figure [T2"l Clearly 
the polynomial fit induces long-range correlations and 
anti-correlations between spectrum bins, and washes out 
the original correlations. We concluded that while we 
could rely on the polynomial fit for the mean values of 
the spectra (since the fit residuals are acceptably small) , 
we could not use the resulting correlation matrices for 
error propagation. 




FIG. 12. Bin-to-bin correlation matrix (25 x 25 bins, 2-8 MeV) 
of converted neutrino spectrum from 235 U, including the ILL 
experimental errors and conversion effects from [2p| ]. 

We therefore choose to propagate the original correla- 
tion matrix on our final binning by interpolating linearly 
between each bin of the original covariance matrix, and 
then computing the resulting correlation matrix. For bins 
which were not present in the original energy range, we 
take a constant error of 20%, fully correlated bin-to-bin 
in order to obtain a conservative upper bound on the er- 
ror. This allows us to use any binning over any energy 
range. Of course the resulting matrix still has rank 25 
over the 2-8 MeV range, but since we do not need its 
inverse this does not limit our ability to use it for error 
propagation. An example of the resulting matrix can be 
found in Figure [TJ] for 235 U. 

Finally, as a cross-check of this method, we computed 
the error on the cross-section per fission for each isotope, 
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FIG. 13. Bin-to-bin correlation matrix (25 x 25 bins, 2- 
8 MeV) of fitted neutrino spectrum from 235 U, following the 
toy Monte-Carlo procedure outlined in the text. While the 
structure of Figure [12] is still present it has been significantly 
washed out and anti-correlation areas have appeared. 




FIG. 14. Bin-to-bin correlation matrix (1000 x 1000 bins, 
1.75-10 MeV) of the fitted neutrino spectrum from 235 U, ob- 
tained after linear interpolation of the original matrix, and 
extrapolating to a fully correlated 20% error outside the orig- 
inal binning. Clearly the experimental correlations from the 
ILL measurements are kept with this technique. 



using our new error propagation method with interpo- 
lated ILL and conversion correlation matrices (column 2 
of table U reproduced in table |V] below), and directly 
folding the converted ILL spectra with the inverse beta 
decay cross-section. Clearly, the results are extremely 
close. The effect of the error on the extrapolated part of 
the spectrum is found to be negligible. We also verified 
that the resulting error bars do not significantly depend 
on the chosen binning. These calculations validate the 
use of the exponential-polynomial fit along with the in- 
terpolation technique for error propagation on the final 



Isotope 


This work 
1000 bins, 1.8-10 MeV 


ILL spectrEt 
25 bins, 2-8 MeV 


prcd 

°/, 235 u 

prcd 
°/, 2 39Pu 

prcd 
f fi 241p u 


6.61±2.11% 
4.34±2.45% 
5.97±2.15% 


6.61±2.13% 
4.33±2.46% 
6.02±2.16% 



TABLE V. Comparison between cross-sections per fission ob- 
tained from fitted spectra (central column), and from a direct 
convolution of the converted ILL spectra (right column). The 
error bars are 'exact' in the case of the direct convolution, 
while they are approximate for the fitted spectrum, obtained 
from interpolation of the original matrices. 



spectrum, and allow us to use the binning of our choice 
without affecting the errors. 

For 238 U, since the spectrum is obtained from ab initio 
calculations, each bin is given an error ranging from 10% 
at low energy to 20% at high energy. However, if these 
errors are uncorrelated, the uncertainty on the overall 
neutrino rate and hence the cross-section per fission for 
238 U is artificially low. To avoid this, we consider that 
there is a fully correlated 10% error on bins from 2 to 
5.5 MeV, another fully correlated 15% error on the block 
from 5.5 to 6.75 MeV and another 20% fully correlated 
error on the last bins of the spectrum. With this con- 
servative prescription our error on the 238 U cross-section 
per fission is 8.2%. 

Finally we derive the total error, crs tot , on the antineu- 
trino rate, for a given core composition. We account for 
the uncertainties on the cross section per fission per iso- 
tope, and the uncertainties on the averaged fraction of fis- 
sion per isotopes. As an example we can consider the case 
ofBugey-4. The fuel composition is given in Ref. Q. The 
new values of the cross sections per fission per isotope 
as well as their uncertainties are given in Table |U The 
fk coefficients are taken as correlated, following a typical 
fuel evolution curve, and such that J2k A-=l±0.6% to ac- 
count for the error on the thermal power of the nuclear 
core. In a first case we consider a 10% relative uncer- 
tainty on the fuel composition. We obtain a final error 
on the expected number of events of 2.7%. Reducing the 
error on the fuel composition to 3.5% results in a final 
error of 2.4%. 



Appendix B: Gallex/Sage, and MiniBooNE 
reanalysis 

In this appendix we briefly provide details on our re- 
analysis of published data from Gallex 51 Cr source data 
SAGE 51Cr and 37 Ar source data [U, and 
MiniBooNE neutrino data [59| . 

The Gallex collaboration performed two measurements 
with two 51 Cr sources in |62j |. obtaining the ratio of the 
measured and expected event rates. These data were re- 
cently reanalyzed [63[, providing updated values of the 
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FIG. 15. (left) Reanalysis of the gallium calibration data of Gallex and Sage. A\ 2 projections as a function of sin 2 (20), without 
correlations (dashed) and with correlations as explained in the text. The dashed contour is very close to that of [56|. validating 
this approach, (right) A\ 2 valleys of the sin 2 (2fl ncw ) parameter of the reactor antineutrino experimental results compared with 
the gallium results. Using the prescription of Ref. [771 ] the parameter goodness-of-fit is 27.0%. 



errors and the ratios. The SAGE collaboration per- 
formed a similar measurement with a 51 Cr source, and 
more recently with a 37 Ar source [6fj| . [6l| . An analysis of 
these data in terms of neutrino oscillation was performed 
in |1[5|. The data values are 0.95 ±0.11, 0.81 ±0.11, 
0.95±0.12 and 0.79±0.09. Performing a neutrino oscilla- 
tion search with these values yields contours very close to 
those of [56| . However we decided to include possible cor- 
relations between these four measurements in this present 
work, that were not previously taken into account. As 
the two Gallex measurements used the same experimen- 
tal technique, we decided to fully correlate their system- 
atic errors, which we understand to be 5.6% and 7.4% 
respectively. The statistical errors remain of course un- 
corrected. For SAGE we followed a similar procedure, 
with systematics of 5.7% and 7.0% according to [6fjl[6l1|. 
Our fractional covariance matrix is 



/ 1.31 
0.41 


V 



0.41 








1.55 











1.53 


0.40 





0.40 


1.30 



10" 



(Bl) 



As expected we obtained slightly more conservative 
contours than those of [56j with this method. 

It turns out that the sin 2 (2# ncw ) best fits ob- 
tained from reactor antineutrino rate results only 
(sin 2 (2# noW! bcstfit) ~ 0.10) and the gallium results only 
(sin 2 (20 neW! bcstfit) ~ 0.27) are slightly different, though 
with large uncertainties. We quantified this effect by us- 
ing the prescription of Ref. [77| . The corresponding Ax 2 
valleys are displayed on the right panel of Figure [T5] Af- 
ter marginalizing over |Am 2 cw (since fitted at the same 
value in both data sets) the parameter goodness-of-fit 
between the two data sets on sin 2 (26> ncw ) is 27.0%, in- 
dicating no significant tension between neutrino and an- 
tineutrino anomalies. 
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FIG. 16. Analysis of the MiniBOONE data following a pro- 
cedure very similar to that of [56|]. Our results are in good 
agreement. 



We also analyzed the MiniBooNE neutrino data 
from (59| following the method outlined in 56] . The idea 
is to fit both electron-like and muon-like data at the same 
time. However, instead of searching for za, — > v e appear- 
ance, we search for v e disappearance, while allowing the 
global normalization of all the samples (e-like and mu- 
like) to fluctuate. 

The excess of e-like data in the low energy part of the 
spectrum is well-fitted by the combined increase in over- 
all normalization and by the disappearance of v e induced 
e-like events. We used a simplified version of the analy- 
sis in reference [56], in that we used a covariance matrix 
independent from the normalization fitting parameter /. 
We obtain very similar results, although our best fit point 
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is in a slightly different location, at sin 2 (2#) = 0.23, We have only included this analysis of MiniBooNE's 
Am 2 = 1.88 eV 2 and / = 1.083, as can be seen in Fig- data for completeness. However in our final statistical 
ureHH significance it has very little impact. 
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